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. . , The influence of multiboson effects on pion multiplicities, single-pion spectra and two-pion cor- 

0^ ' relation functions is discussed in terms of an analytically solvable model. The applicability of its 

basic factorization assumption is clarifled. An approximate scaling of the basic observables with 
, the phase space density is demonstrated in the low density (gas) limit. This scaling and also its 

' ' violation at high densities due to the condensate formation is described by approximate analytical 

^ I formulae which allow, in principle, for the identification of the multiboson effects among others. 

O . For moderate densities indicated by the experimental data, a fast saturation of multiboson effects 

' with the number of contributing cumulants is obtained, allowing for the account of these effects in 

realistic transport code simulations. At high densities, the spectra are mainly determined by the 
universal condensate term and the initially narrow Poisson multiplicity distribution approaches a 
' wide Bose-Einstein one. As a result, the intercepts of the inclusive and flxed-n correlation func- 

I , tions (properly normalized to 1 at large relative momenta) approach 2 and 1, respectively and their 

^ ' widths logarithmically increase with the increasing phase space density. It is shown that the neglect 

ijr^ ' of energy-momentum constraints in the model is justified except near a multipion threshold, where 

1/^ , these constraints practically exclude the possibility of a very cold condensate production. It is ar- 

' gued that spectacular multiboson effects are likely to be observed only in the rare events containing 

sufficiently high density (speckle) fluctuations. 
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^ ■ I. INTRODUCTION 

^ : 
^' 

, In future heavy ion experiments at RHIC and LHC one expects to obtain thousands of pions per a unit rapidity 
■ interval. Since the pions are bosons there can be multiboson effects enhancing the production of pions with low relative 
. ^ momenta thus increasing the pion multiplicities, softening their spectra and modifying the correlation functions. 

Though the present data does not point to any spectacular multiboson effects, one can hope to observe new interesting 
I phenomena like boson condensation or speckles in some rare events or in eventually overpopulated kinematic regions 
with the pion density in the 6-dimensional phase space, / = {2n)'^d^n/ dPpdPx, of the order of unity (see, e.g., 

B)- . . . 

In the low-density limit (/ <^ 1), the mean phase space density at a given momentum p can be estimated as 
the mean number of pions interfering with a pion of momentum p (rapidity y and transverse momentum pt) and 
building the Bose-Einstein (BE) enhancement in the two-pion correlation function (/)p ~ 'K^^^N{p)/V, where 

N{p) = (Pn/d^p and V = r^ryrz is the interference volume defined in terms of the outward (r^), sideward (ry) and 
longitudinal (r^) interferometry radii. Typically (/)p ~ 0.1 for mid-rapidities and pt ^ (pt) 0|- The data are also 
consistent with the phase space density of pions near the local thermal equilibrium |^ . 

At AGS and SPS energies the interference volume V seems to scale with dn/dy (see, e.g., [p^|jTl| ]) pointing to the 
freeze-out of the pions at a constant phase space density.^ If this trend will survive then there will be no spectacular 
multiboson effects in the ordinary events at RHIC or even at LHC. In such a situation the standard two-particle 
interferometry technique could be used to measure the space-time intervals between the production points also in the 



^Similar effect was observed also for protons produced in hadron- and electron- nucleus interactions p^ . 
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future collider experiments. The corresponding interferometry radii for lead-lead collisions at LHC would be however 
rather large - about 20 fm. 

The multiboson effects can show up however in certain classes of events. An example is a rapidly expanding system 
with the entropy much smaller than in the case of total equilibrium. Then a strong transverse flow can lead to rather 
dense gas of soft pions in the central part of the hydrodynamic tube at the final expansion stage (see, e.g., |p^). 
Another reason can be the formation of quark-gluon plasma or mixed phase. Due to large gradients of temperature 
or velocity the hydrodynamic layer near the boundary with vacuum can decay at a large phase space density and lead 
to pion speckles even at moderate transverse momenta [Q. 

The dramatic difference in behavior of Boltzman-like gases and dense multiboson systems can lead to serious 
problems for transport models like RQMD, VENUS, etc., ignoring actually the statistical properties of the particles 
both in intermediate and final states. In these models the most intensive particle production happens at relatively 
early evolution (expansion) stage when rather large pion phase space densities can be achieved at RHIC or LHC 
energies. 

Generally, the account of the multiboson effects is extremely difficult task. Even on the neglect of particle interaction 
in the final state the requirement of the BE symmctrization leads to severe numerical problems which increase 
factorially with the number of produced bosons [^[^. In such a situation, it is important that there exists a simple 
analytically solvable model Q allowing for a study of the characteristic features of the multiboson systems under 
various conditions including those near the Bose condensation. In this paper we use this model to demonstrate 
the infiucncc of the multiboson effects on pion multiplicities, spectra and two-pion correlation functions. Besides 
the original papers , similar studies can be found also in p5|-[l7|] . Particularly, some of the new aspects of the 
multiboson effects, like the scaling behavior with the phase space density or the behavior of the (semi-) inclusive 
correlation functions near the condensation limit were studied in our unpublished paper |T5| . The present work 
represents an elaborated version of the latter. 

In Section || we introduce the space-time description of particle production in terms of Wigner-like densities and 
discuss their physical meaning and the conditions of their factorization in the model of classical one-particle sources. 



The multiboson formalism in the factorizable case is set forth in Section III. Using this formalism and the simple 
Gaussian ansatz for the emission function, we present in Section IV the analytical solutions (partly in terms of the 
recurrence relations) for the multiplicity distribution, single-boson spectra and two-boson correlation functions. In 
Section we compare the results of numerical calculations with the analytical approximations accounting for the 
approximate scaling behavior in the low density (gas) limit as well as for the condensate formation at high densities. 
The results are discussed and summarized in Section Vl. 



II. SPACE- TIME PICTURE OF PARTICLE PRODUCTION 



A. Wigner-like density 

Let us first consider a process in which, besides others, just n non-identical particles of given types are produced 
with the 4- momenta Pi = {Ei,pi} and Lorentz factors ji — Ei/nii (to simplify the notation, we assume that particles 
are spin-less). The inclusive differential cross section of this process is described by the invariant production amplitude 
Tn{pi, . . . ,Pn;a): 

^1 ■ • •T"l3 ~Pi = XI I Tn{Pl, ■ ■ ■ ,Pn;a) pEE 71 • • • 7„Cr„P„ (pi , . . . ,P„), (1) 

where the sum over the quantum numbers a, describing the rest of the produced system, contains also an integration 
over the momenta of the other produced particles with the energy- momentum conservation taken into account. The 
non-invariant production probability P„(pi, . . . ,p„) is normalized to unity. 

If the particles are identical spin-less bosons, then the production amplitude has to satisfy the requirement of Bose 
symmetry. Formally, this can be achieved by the substitution of the non- symmetrized amplitude T'„(pi, . . . ,p„;a) 
(corresponding to the "switched off' effect of quantum statistics) by a properly symmetrized one: 

Tn{pi,...,Pn;a) {niy^/^^Tn{pai, ■ ■ ■ ,Par,;a), (2) 
where the sum is over all n\ permutations a of the sequence {1, 2, . . . n}. 
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In the following we will neglect particle interaction in the final state.^ Then the non-symmetrized amplitude 
Tn{pi, ■ ■ ■ ,Pn] oi) = T„(V; a) is related to the amplitude in the space-time representation Tn{xi, . . . , x„; a) = Tn{X; a) 
(describing production of the particles in the space-time points with the 4-coordinates A" = {xi, . . . , a;„}) by the usual 
Fourier transform: 

Tn{V;a)= j dX(^'^^Tn{X;a). (3) 

Inserting @ into and introducing the space-time density matrix pn{X\X') = ^Tn{X\Q)T*{X';a) and its 
partial Fourier transform (emission function) - an analogy of the Wigner density |20[| : 

D„ipi,xi; . . .■,pn,Xn) = 71 • • •7nO'ri-D«bi,ii; • ■ • ;p„,x„) = / dSe'-'^^ pn{X + }^5\X - \5), (4) 
where X = ^{X + X'), S — X — X' , we can rewrite the production cross section (0) in the absence of BE effects as: 

„ ^ '^"o = cr„P„(pi, . . . ,p„) = 0-„ / d'^Xi- ■ ■d'^XnDn{pi,Xi; . . .■,Pn,Xn)- (5) 

Since Z?„ is a real (though not positively defined) function normalized to unity, in accordance with Eq. (|^) it can be 
considered as an approximation to the emission probability of the particles with given 4-momenta pi in the average 
space-time points Xi = ^{xi + x'^). 

The insertion of the symmetrized amplitude (||) into the cross section formula (0) leads to the substitution of the 
probability P„ by a BE modified one P^. For example, in case of a two-boson production process, instead of Eq. (|^) 
we have: 

^'2(Pi,P2) = J d*xid'^X2[D2ipi,xi;p2,X2) + D2{p,xi;p,X2)cos{qx)], (6) 

where p — \{pi + P2), 9 — Pi — P2 and x = xi ~ X2. Clearly, for n > 1 the probability is no more normalized to 
unity. The integral over this probability yields the BE weight ujn of an n-boson event produced in the absence of the 
effect of quantum statistics (o^o = o^i = 1, > 1 for n > 1): 

y" rf3pi---rf3p„P,^(pi,...,p„)/ j d^Pi---d^P„P„(pi,...,P„). (7) 

We will also use the differential BE weights uj'if'{'Pi, ■ • ■ ,Pfe) defined as in Eq. (^ except for a skipped integration 
over the momenta of particles 1, 2, . . . , fc. 



The multiboson problem greatly simplifies (see Section III) in the factorizable case when the n-particle emission 



function and, as a consequence, the non-symmetrized production probability can be written as products of the single- 
particle ones: 

Dn{pi,Xi; . . . ■,Pn,X„) = D{pi,Xi) ■ ■ ■D{pn,Xn), Pn(pi, . . . , p„) P(pi) ' ' ' P{Pn)- (8) 

Consequently, the BE weights are expressed through a universal function Fij (see, e.g., p8[): 



4"nPl,P2,...Pn)=En^-.' (9) 

where 

F,, = J d''xD{p,,,x)-exp{iq,,x)/ [P(p,)P(Pj)]'^', (10) 

Pij ~ \^Pi '^Pj) lij — Pi ^ Pj- The sum in Eq. (|^) is over nl possible permutations a of the sequence {1,2, . . . n}. 
For example, in the two-boson case, we have 

'^2'^(Pl,P2) =^^11^^22+^12^21 = l+i^f^(pi,P2), (11) 



(2) h -H 

where K^ ' = P12P21 is so called differential cumulant (see Section |I|) . 



^This is more or less valid assumption for neutral pions but not for the charged ones. For the treatment of multiboson effects 



in the case of interacting pions see, e.g., Jlit . 
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B. One-particle sources 



To clarify the physical meaning of the emission function and the factorization assumption (H), let us follow Kopylov 
and Podgoretsky (see, e.g., 0) and assume that particles 1,2,... are emitted by one-particle sources A,B,... which 
are considered as classical so they can be treated by parameters and not by amplitudes. Thus the 4-coordinates 
of the source centers xa,xb, ■ ■ . and other source characteristics in the model can be considered as a part of the 
quantum numbers a = {xa, xb, ■ ■ ■ a'}. It was pointed out by Kopylov and Podgoretsky that the BE effect is mainly 
determined by the phase factor expiipiXA + ip2XB + • • •) contained in the amplitude T'„(pi, . . . ,p„; a). 

Let us first consider the production of only one boson. Assuming the translation invariance of the decay amplitudes 
u, we can write the single-boson amplitude in the 4-coordinate representation as: 

T(xi; a) = t{xi — xa] ol) = u{xi — xa] oi)v(oi), (12) 

where a = {xa, a'} and a' = {A, . . .}. Inserting (|l2|) into (^) and introducing the Fourier transform 

t(p; a)= ( d^^e'PH{£,] a) = u{p] a')v{a), (13) 



we obtain the Kopylov-Podgoretsky ansatz: 

T{p] a) = e'P'^^tip] a) = e'P'^^uip] a')v{a). (14) 
For the production probability we have 

P^(p) = P(p) = I d^XAWip, xa) - J d'xA I t{p] {xA, a'}) |2, (15) 

a' 

where we have introduced the (true) emission probability W{p,xa). Similarly, if two bosons are produced, the 
production probability takes on the form: 

Pi{Pi,P2)= J d-^XAd'^XBj2{\tiPi^P2]ct) P -l-i?e[t(pi,p2;a)i*(p2,Pi;a)e*«^]} . (ig) 

a' 

Note that here x = xa — xb and, as usual, q — pi — P2- 

The emission function D can be expressed through the Kopylov-Podgoretsky amplitudes t{p] a) continued off mass- 
shell. Using the inverted Fourier transform in Eq. (|l^), we get from Eq. (Q) |18|] : 

D{p, x,)^j^J:J d^XAd^Ke'^^^^-^^') ■ t{p + \k] {xa, a'})i*(p - \k] {xa, a'}) 

= S / d^x^d^^e-f-^-^i) • u(j, -f \k] a')u*{p - \k] a')\v{a)\\ ^^^^ 

a' 

It is clear from Eq. ( p^ ) that the momentum dependence of the amplitude t[p] a) is determined by the space- 
time extent of the one-particle source. For example, assuming that the source emits a particle independently of the 
quantum numbers a' except for the source type A: u(^] a') = 'u(^; A) and that the distribution of the emission points 
in the source rest frame is given by a simple Gaussian with the width parameters ta and ta characterizing the proper 
space-time sizes of the source A, we obtain in case of a source at rest 

t{p] a) ~ exp(-ir^p2 - ^Tlpl)v{a). (18) 

The probability |i;(q;)|^ describes the production of particle sources and depends on the 4-coordinates xa of the 
source centers. In the following we will take it also in a simple Gaussian form with the width parameters Fq and tq 
characterizing the space-time extent of the source production region. 



^Note that Eq. (|l|) is valid also off mass-shell when po E. For a source moving with a non-relativistic velocity /3^, the 
substitution p ^ p — PA has to be done in Eq. (|l^), where pA = is a mean 3-momentum of the particle emitted by a 

source A. 
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Comparing Eqs. (|1^) and (|17|), it can be seen that the emission function D is more spread in space and time than 
the emission probabihty W. In particular, the Gaussian parameterizations of A) and |u(a)|^ yield: 

2 2 

Wip, x) ~ eM-r\p' - tIpD e^-^ " (19) 

and 

2 2 

D{p,x) ~ exp(-r^p2 ~ rlpg)exp(- ^ - a'^" J . (20) 

Clearly, the factorized form (^) for the emission function is recovered in case of independent sources (i.e. sources 
having no quantum numbers in common) assuming a unique mechanism of their production. Generally, the latter 
condition may not be fulfilled, e.g. in case of heavy ion collisions without selection of the impact parameter. Then, 
even for independent sources, Eq. (^) will be substituted by a weighted sum of factorized terms corresponding to 
different single-particle emission functions. 

The single-particle density generally contains contributions from the sources of different type {e.g. different res- 
onances). It is interesting to note that in case of only one source type {i.e. universal source parameters) and on 
condition of sufficiently slow relative motion of the sources contributing to low-|g| pairs {e.g., due to limited source de- 
cay momentum) , the BE correlation effect in the two-boson case is solely determined by the characteristic space-time 
distance between the source centers: 

P|(pi,P2) =P(pi)P(p2)[l + exp(-?^q2-f2q2)]_ (21) 

This result follows from Eq. (^) or, in terms of the Kopylov-Podgoretsky emission amplitudes, immediately from Eq. 
(H). If the sources of interfering bosons move with non-relativistic velocities = p^/m and the distribution of 
their characteristics is given by a product of exp(— p^/2Aq) and the same Gaussian as before {i.e. the residual slow 
relative motion decouples from other source characteristics), we still arrive at Eqs. (^9|), (pO|), up to a substitution 
f\v^ ~* f^P^ /^{"''A^oY' + 1] corresponding to a widening of the momentum distribution due to the source relative 
motion. Note that the factorized form of the multiparticle density is not destroyed by this motion. The latter however 
influences the two-boson correlation function which now becomes sensitive to the source size ta even in case of one 
source type, Eq. ( ^l| ) being modified by the substitution 7q ^ ?o + ''a/I^ + (^aAq)"^]. 

The space-time extent of one-particle sources can be usually considered much smaller than the characteristic space- 
time distance between their centers (r^ <C rg, ta ^ tq). The 4- momentum dependence of one-particle amplitudes 
is then negligible when varying the particle 4- momenta by the amount ~ ^o'^i'^o'^ characteristic for the interference 
effect. On such a smoothness condition, there is practically no difference between the emission function D and the 
emission probability W and both Eqs. ^ and (|6|) yield the well known result of Kopylov-Podgoretsky for the 
production probability of two identical bosons: 

^'2(Pl,P2) = ^^(Pl,P2)[l + {C0s{q{x^ - X^)))]. (22) 



III. MULTIBOSON FORMALISM IN FACTORIZABLE CASE 



The multiboson effects can be practically treated provided that we can neglect particle interaction in the final 
state and assume independent emission of non- interfering particles (a valid assumption for heavy ion collisions), 
supplemented by the requirement of a universal single-particle emission function D{p, x) for the detected class of 
events. We can then use Eq. (H) expressing the n-particle emission function as a product of the single-particle ones. 
Then, similar to refs. |^,|| it is convenient to define the functions 

Gi(pi,P2) = J d'^xD{^{pi +P2),x) ■ exp{i{pi -p2)x), 

G„(pi,P2)= j d^k2...d^k„Gi(pi,k2)...Gi(k„,p2) 

EE j d3k2G„_i(pi,k2)Gi(k2,P2), 

5„ = y" d3pG„(p,p). (23) 
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The function Gi at equal momenta is just the initial (not affected by the multiboson effects) single-boson spectrum 
normalized to unity: 



P(p) = Gi(p,p), J d^pPip) = <?! = 1. 
The related quantities are so called cumulants 

n-l 

Ki'HpuP2) = (n-2)!^G.(pi,p2)G„_.(p2,pi)/ [P(pi)P(p2)], 



(24) 



i=l 



i^(i)(p)-(n-l)!G„(p,p)/F(p), 
Kn = [n- ly.gn- 



(25) 



It can be shown that the BE weight of an event with n identical spin-zero bosons is determined through the cumulants 
Kj by the recurrence relation : 



W„ = Gq ^KiUJn-l + G" '^K20Jn- 



n— 1 j 



(26) 



with ujf) = uJi = 1; C"~^ = [n — l)!/[i!(n — 1 — i)!] are the usual combinatorial numbers. For example, ^2 = 1 -I- K2 
and W3 = 1 + 3A'2 + A'3. One can check that w„ — n\ provided that all the elementary one-particle sources are 
situated at one and the same space-time point so that all the single-boson states are identical and ifj+i = In the 
other extreme case of a large phase space volume and n?'K2 <C 1, we can neglect the contribution of the higher order 
cumulants except for the first power of K2 and write: 

UJn = l + C^K2. (27) 

Given the initial multiplicity distribution u)(n), the BE affected one is easily calculated using the BE weights aj„: 

00 

w{n) = ujnw{n) I ^ujjw{j). (28) 

j=o 

Particularly, assuming the initial Poissonian distribution with the mean multiplicity rj: w{n) — e^'^rf^ /n\, we get: 



(n) =cj„— -/ > t^j-T. 



(29) 



Similarly, the BE affected single- and two-boson spectra, respectively normalized to n and n[n— 1), can be written as 

7V(i) (p) ^ iV„(p) - n^i) (p) • P(p) / a;„ (30) 

and 

iV(2Hpi,P2) =n(n-l)42)(p^,p2)P(p,)P(p2)/c.„, (31) 



where the differential BE weights ti;i^-'(p) and Li;^''''(pi,p2) are expressed through the differential cumulants K)^^ 
andxi')(pi,P2) lH: 



,(2) 



(i; 



n-2 



/n— 1 
d'p'J^\p,p')Pip') = ^G;-iif«i(p)a;„_i,„ 



OJ, 



^ G/XW (pi)A:ji\+i(P2) + Kj%{PuP2) 



.1=0 



(32) 



*This situation is similar (flat correlation function) though different from the case of the emission of so called coherent bosons 
for which there is no enhancement factor. In fact, when the one-particle sources become closer and closer, so that their distances 
are less than the wave length of the emitted bosons, they can no more be considered as independent ones and a multiparticle 
source of non-interfering bosons has to be introduced [Ell. To quantify the transition to the non- interfering bosons a concept 
of the coherence length can be used M . 
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The differential weight cun (pi, P2) can be considered as a two-particle correlation function measuring the BE effect 

— (2) 

on the initial uncorrelated two-particle spectrum Nn (pi,p2) — ^(^ ~ l)^(Pi)^(P2)i with the normalization 



/■ 



d''pid''p2iol^\pi,P2)P{pl)P{p2) = iOn. (33) 

Usually the correlation function is normalized to unity at a large |q|. Such a normalization is approximately satisfied 
for the correlation function defined as: 

R„{p,,P2) ^ N^^'Hpi,P2) / Ni^\pi,P2) = L^l?\pi,P2) / co^. (34) 
In practice, the two-particle correlation function is defined through the observable spectra as: 

i?n(Pl, P2) = CrX?\pi,p2) I {N\l\p{)N\l\p2)\. (35) 
Similarly, the (semi-)inclusive correlation function is defined as 

i?(pi,P2)=ciV(2)(pi,p2)/ [iv(i)(pi)iv(i)(p2)], (36) 

where 

A^W(p)=^«;(n)7VW(p), iV^^) (p^, p^) = ^ ^(n)7V(2) (p^, p^) (37) 

n n 

are the corresponding (semi-)inclusive single- and two-particle spectra, w{n) is the normalized multiplicity distribution 
accounting for the BE effect according to Eq. (p8|). Later on, using an analytical Gaussian model for the emission 
function, we show that the normalization constant c„ can be expressed through the BE weights as: 

c„ = nuj\_^ I {(n - l)tj„cj„_2] (38) 

and that c — \ for the inclusive correlation function provided a Poissonian multiplicity distribution of the initially 
uncorrelated bosons. 

As one can see from formulae (pO[)-(^), the multiboson correlations lead to distortions of the initial single- and 
two-particle distributions. Such distortions are small in the case of interference of only two or three identical particles. 
However, they can become essential for the events with a large number of identical bosons due to factorially increasing 
number of the correction terms Q (see also 0] and |^). For the processes characterized by a high (> 0.1) phase space 
density of the identical bosons at the freeze-out time the multiboson effects can no more be considered as a correction 

To account for the multiboson symmetrization effect in the event simulators, a phase space weighting procedure 
was used with weights in the form of a normalized square of the sum of n\ plane waves This procedure however 

appears not practical for a large n due to the factorially large number of the terms to be computed to calculate the 
weight and, due to large weight fluctuations. These fluctuations can be substantially reduced by weighting only in 
the momentum space. The corresponding BE weights are given in Eq. (|p). They are expressed through the universal 
function ( |lO| ) which is simply related with the function G\. 

F,, = Gi(p.,p,) / [P{p,)P{p,)Y/\ (39) 

On the condition of sufficient smoothness of the single-particle spectra, we can put 

Fij ^ j d'^xD{pij,x) ■ exp{iqijx) / J d'^xD{pij,x), (40) 

where pij = ^{pi + Pj) and = pi — pj. This function can then be calculated as suggested in |l^: 

Fij = {cxp{iqijXk))p^^ , (41) 

where the averaging is done over all simulated phase space points {p^, Xk} such that pj, is close to a given 3-momentum 
Pij. However, there is still the problem with factorially large number of the terms required to calculate the weight 
according to Eq. (||). 
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Fortunately, when calculating only single- or two-particle distributions according to Eqs. I pUi ) or (^), this number 
is strongly reduced (eaten by the combinatorial numbers C™ in Eqs. (^^). We should however perform integration 

over momenta of one or more particles to determine the integrated cumulants i^i^''(pi, P2), Kn\p) and Kn- 

The numerical averaging of the cumulants of all orders is a difficult task. In case of large multiplicities of identical 
bosons {n > 20) this is practically possible in the models with a symmetric emission function (allowing to use a 
special Monte Carlo technique) [|j or with a simple analytical parameterization of this function [gH]. For example, 
in ref. Q the corrections to multiplicity distributions, single-particle spectra and two-particle correlation functions 
were calculated using the relativistic Bjorken model [p3| for the emission function. To compute cumulants up to tenth 
order, the integration was performed analytically over the space-time coordinates and numerically over the momenta. 

Generally, for realistic models used to predict particle production in ultra-relativistic heavy-ion collisions, the 
numerical limitations allow to determine only a few lowest order cumulants (up to about the fourth order) [|l8| . 
Fortunately, since the interferometry measurements point out to a moderate pion freeze-out phase space density of 
^0.1, the lowest order cumulant approximation appears to be reasonable for typical events in present and likely also 
in future heavy-ion experiments (see Section ^). At the same time, even in the absence of strong multiboson effects, 
their account can still be important for realistic simulations of heavy ion collisions fl^ . 



IV. ANALYTICAL MODEL 



To study the multiboson effects in a dense pion gas, we use a simple model assuming independent particle emission 
(see Eq. (^)) with the Gaussian ansatz for the single-boson emission function D{p,x) 

^(^'^)-(2^^^P(-2Z^-^)^W. (42) 

Note that this ansatz corresponds to the independent one-particle sources of Kopylov and Podgoretsky, all of the 
same type, characterized by a universal size of ~ (2A)~-'^, with the centers distributed according to a Gaussian 
of a dispersion = — (2A)~^ (see Eq. (pO|)). Then, in the low density limit but regardless of the validity 
of the smoothness condition ro ^ (2A)~^ (see, however, the footnote after Eq. ( |2^ ) concerning the independence 
assumption), the correlation function of two non-interacting identical particles measures the dispersion of the relative 
4-coordinates x of the centers of the one-particle sources as the inverse width squared of the correlation effect seen in 
the relative momenta q = Pi — P2 For spin-0 bosons 

^(Pi, P2) = 1 + (cos(qx)) = 1 + cxp(-?^q2). (43) 

In this model the initial boson phase space density (not affected by the BE effect) is given by 

~ n p"^ 

^"(P'")=(W^"P(~2A^-2^)- 

The mean densities at a fixed boson momentum p and averaged over all phase space are 

(/„)p - / d^^iUr I J rf^x/„ = ^^=^ exp(-^) (45) 

and 

ifl) = Jd^xd^pifnf I J d^xd^p/„ = n/ (2roA)^ (46) 

respectively. Similarly, the initial inclusive densities /(p,x), (/)p and (/) are given by Eqs. (Q)-(^) with the 
multiplicity n substituted by the initial mean multiplicity. 

It is worth noting an approximate equality (see also a model independent prove in Q) between the mean phase 
space density in the low density limit 

(/>P = exp(— ^) = ^iV(p) (47) 

(V2roA)3 ^ 2A2^ ^ ' 

and the mean number of pions building the BE enhancement in the two-pion correlation function 
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/ d^ci[RiPi,P2) - l]N{p,)N{p2)/N{p) = ^N{p) « (/)p, (48) 



r3/2 

'o 

Pi 2 = P ± q/2. This equality is valid up to relative corrections 0((roA) ~^) and 0((/)p), the latter representing an 
impact of the BE correlations on the single-boson spectrum (see Section VB ). 

It is important that the Gaussian ansatz in Eq. ( ^ ) allows to express the functions G,i(pi,p2) and the integrals 
gn (see Eqs. (23)) in simple analytical forms ||l5[] : 

G„(pi, P2) - {2TTA'An)-^/' exp(-&+(pi + p2f - 6-(pi - P2f), 

g„ = (8A2A„6+)-3/2, (49) 

where A„, 6^ and b~ are given by the recurrence relations: 

An = 2A2A>-i(6+_i + +b+ + b^), 

b-=b+b-/b+, (50) 

with Ai = 1, b^ = 1/(8A^) and b^ = ^o/^- The recurrence relations of this type allow for the analytical solution 
prf . In our case it reads as: 

6+ = 5+6-1(1 -p")/(l + p"), ^„ = - P'"), 5n=r"(l-p")-', (51) 
where e"^ = 2roA, p = (1 — e)/(l + e) and the parameter 

/3 = (roA + 1/2)3 (52) 

can be considered as a characteristic phase space volume. 

For example, for n = 2 and 3, we have: A2 = (1 + e^^)/2, A3 = (1 + Se-'^^/W, b+ = 26+/(l + e^), b+ = 
36^(1 + e^/3)/(l + 3e^), 512 = e^, 53 = [4e^/(3 + e^)]^- Recah that the cumulants related to 32 and 53 are (see Eq. 
(||) K2 = 92 and K3 = 2g3. 

It follows from the recurrence relations (|o|) or their analytical solutions ( |5l| ) that the slope parameters 6+ and b~ 
approach each other with increasing n. In the large-n (n > tqA) limit we then have [p^ : 

5+ ^ 5- ^ ro/(4A), A. - (2roA), 5„ ^ /?" " (53) 

and 

G„(Pi,P2) - /3-" (^)'^'exp (-^(4p^ + q')) ■ (54) 

In very large-n (n > e(3) limit, using the large-n behavior of the parameters g„, we can get from the recurrence 
relation (26) the following behavior of the BE weight (isj: 

c^„^c(/3)n!//3", (55) 

where c(/?) is a function factorially increasing with f3, c(l) = ij^ 

It is worth noting that the large-n limits become equalities a.t f3 — 1 (rgA = 1/2) when gn = An = 1, A'„ — {n— 1)!, 
LUn = n! and 6,7 = 6+ = ro/(4A). Recall that (3—1 corresponds to the minimal possible phase space volume when 
all the particle emitters are situated at one and the same space-time point so that the size (2A)~-'^ of the elementary 
source determines not only the width of the single-particle spectrum but also the characteristic distance between the 
production points (see however the footnote after Eq. (|2^) and also ref. |18| for a more detailed discussion). In such 
a case ro = and the correlation function equals 2 for any value of q. 

In the low-n (n < tqA) limit, i.e. in the case of a large phase space volume, it follows from Eqs. (|5^) or ( |5l| ) 
that the slope parameter 6+ increases linearly with n up to the corrections 0((2roA)-^) and that, at n > 2roA, this 
increase saturates at ro/(4A). Similar behavior shows the parameter A„//?^"/'^. Thus, at n <C 2roA, we have: 



good approximation is c(/3) = d(/3) = ai + a2/3''^ ai = 0.617, 02 = 0.621 and 03 = 0.788. 
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6+ = = n6+ ^„=n(roA)2(»-i), 5„ = n-3(roA)-3("-i) . (56) 



Comparing the low-n approximations (pq ) for the parameters b^, An and g„ with the large-n ones in Eqs. (p3|), we 
can see that they tail each other at n = 2roA, roA/2 and roA respectively. Correspondingly, the low-ri approximation 
for the G„-function 

G„(pi, P2) = ro3(roA)-3"(2™)-3/2 cxp(-p2n/2A2 - q2roV2n) (57) 

tails with the large-n one in Eq. ( |5^ ) at n = rit, 1/2 < nt/(roA) < 2. 

Consider now the correlation function i?„ defined in Eq. (35). To determine the normalization constant c„, it is 
convenient to rewrite the single- and two-boson spectra at a fixed multiplicity n as 



iPj - l^j^Q ^J^^ Gj + i(p, pj = Gj + i(p,pJ, 

_ i^,i-2-j7("-2-j)! v^j 



jvr(Pi,P2) = j:;zo " ELo[g'+i(pi^pi)g.-'+i(P2.P2) + (gg) 

G; + l(pi,P2)Gj-; + i(p2,Pl)] 



ELo '"'^^fa) ELo[G;+i(Pi,Pi)G'j-;+i(p2,P2) + G';+i(Pi,P2)G'j-i+i(P2,Pi) 



where w{n), defined in Eq. (|29|), coincides with the BE affected multiplicity distribution arising from the Poissonian 
one characterized by the initial mean multiplicity rj and Gi(pi, P2) = ?7'Gi(pi, P2). Noting further that approaches 
the hmiting value ro/(4A) from below, while b~ does it from above, we can see from Eq. ( ^ that, at large q, all 
terms in Eqs. ( |5^ ) for iV^"^^(pi_2) and Nn\pi,P2) (Pi,2 = P± q/2) can be neglected except for those containing the 

lowest slope b^. For the normalization constant c„ = limg.^oo[^n^''(Pi)-^n^''(P2) / -^n'^''(Pi,P2)] in Eq- ( ^Sf ) for the 
correlation function we thus get ||l5| : 

c„ = [w{n - 1)]^ / [w{n)w{n - 2)] = nujl_^ / [{n - l)tj„tj„_2]- (59) 

Note that C2 — 2/aj2 G (li2); with the increasing multiplicity c„ decreases and, according to Eq. (|55|), c„ ^ 1 
for n > e/3. For large phase space volumes (when = 1 at small n), the normalization c„ = n/{n — 1) and the 
exclusive correlation function i?„, normalized to 1 at large |q|, becomes close to the usual definition as a ratio of the 
two-particle spectrum to the product of the single-particle ones, both spectra normalized to 1. Generally, the latter 
definition is however not reliable since it leads to the plateau height of w„a;„_2/w^_i > 1 which, in case of a small 
phase space volume f3 and a small n, can be substantially larger than 1. For example, for n = 2 this height is UJ2 and 
can reach a value of 2 if /? 1 . 

Regarding the (semi-)inclusive single- and two-boson spectra, they can be written in a form similar to Eqs. (|5^ ) 
only in the initially Poissonian case: 

n-l 

^^'^ (p) - E E - 1 - j')^^^! (P' p) / E ^(")' 

n J— n 

n-2 3 ^ _ 

A^^'Hpi, P2) = E E - 2 - i) E[^'+i(Pi' Pi)^J-'+i(P2' P2) + 

n j=0 l=Q 



Gi+i(pi,P2)Gj_i+i(p2,pi)] / E^(^*)- (60) 



The normalization constant in Eq. ( p6[ ) for the (semi-)inclusive correlation function is then |15|: 

c = E - 1)]' / E E "^(^ - 2)]- (61) 

n n n 

Clearly, in the completely inclusive case (when the sums include all n from to oo and w{n — j) = 1), we have 
c — 1 and 

oc 

("> = E = 5, 
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3=0 



oo J 



= G(pi,pl)G(p2,P2) + G(pi,P2)G(p2,Pl), 

where g„ = 77" 5„ . For the inclusive correlation function we have 

G(pi,P2)G(p2,Pl) 



^(Pl,P2) = 1 



G(pi,Pl)G(p2,P2) 



(62) 



(63) 



Thus, in the considered case of the initially Poissonian multiplicity distribution, the intercept i?(p, p) = -R(O) = 2 
in agreement with the result generally valid for thermalized systems Note that Eq. (|6^ ) coincides with similar 
expressions in refs. |16 1^ up to a normalization factor (n)^/(n(n — 1)). With the increasing density, the latter 
decreases from 1 to 1/2 and, for dense systems, forces the corresponding correlation hmction to 1. Such a behavior 
was incorrectly interpreted as a coherent effect (see also a discussion in Section VI). 



V. RESULTS 



A. Multiplicity distributions 

We will consider here the multiplicity distribution (p9|) resulting due to the BE effect on the initially Poissonian 



one with the mean multiplicity rj. In accordance with the large-n behavior of the BE weights in Eq. (55), it takes on 
the following limiting form at n > e/3 jist : 

w{n) -> const' -C, C = v/f^- (64) 
The large-n behavior of the multiplicity distribution in Eq. ( |6^ indicates that it approaches the BE one: 

u;BE(n) = z^"/(l + i^r+\ 1^^^/(1-0, (65) 

with the mean multiplicity v.^ This is demonstrated in Figs. 1 and 2. Thus, at tq = 2.1 fm and A = 0.25 GeV/c, 
the BE effect transforms the initial Poissonian multiplicity distribution with 77 30 (dotted curve in Fig. la) to the 
one with much higher mean and dispersion values (solid curve in Fig. la). The exponential tail expected for the BE 
distribution is clearly seen in Fig. lb where the results are presented in logarithmic scale for 77 = 10, A = 0.25 GeV/c 
and ro = 1.5 fm. One may see that Eq. (joj) (dashed line) becomes an excellent approximation for n > 30, which is 
close to the condition n > ef3 = 37.6 for the present choice of parameters. The slope parameter b in the exponential 
fit w{n) = const • exp(— 6n) of this tail at large n should be, according to Eq. (|6^), only a function of the variable ^: 
b = — In(^). Such a scaling is demonstrated in Fig. 2a for various values of 77, A and tq. Note that ^ = 0.95 and 0.72 
for Figs, la and lb, corresponding to 6 = 0.02 and 0.27, respectively]^ 

It should be noted that the experimental data point to a moderate value of the density parameter ^. Thus, taking 
0.2 as an estimate of the inclusive phase space density at p = from AGS and SPS experiments and using Eq. (jij 
we get (see the last Section) ^ « 0.4 — 0.5. 



This is in accordance with the appearance of the Bose-condensate in a dense ideal Bose gas The fluctuations of the 

number of particles in the condensate are very large - they are described by the well-known Einstein formula for identical 
bosons in the same quantum state. The corresponding BE multiplicity distribution in Eq. (^5|) turns to the Reley one for very 
large mean multiplicities. This type of BE condensate should not be mixed up with the multiboson coherent (laser) state in 
which the BE correlations are absent and the multiplicity distribution corresponds to the Poisson law. 

^ At the explosion point ^ = 1 the tail of the multiplicity distribution becomes a constant {b = 0) so that the mean multiplicity 
(n) would go to infinity provided that there are no energy-momentum constraints. Note that the corresponding critical initial 
mean multiplicity rjcr = /3 = (roA -I- 1/2)^ is close but different from that given in Eq. (9) of ref. For the origin of this 



difference see discussion in Section VI 
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Using Eg. ( |62| ) for (n) and tailing the large- and small-n approximations of the integrals gn at rit = vqA (see Eqs. 
(^) and (|5^)), we can approximate the mean multiplicity as 



= (">g + (">c, (66) 

where ^ = r//(7-oA)^ > ^; the density parameters ^ and ^ coincide at roA ^ 1. At large phase space volumes, 
(roA)'^ '3> 1, the two terms in Eq. ( |6^ ) can be considered as contributions of the BE gas and BE condensate 
respectively. It can be seen that the condensate dominates on condition (n) > /?. 

Note that in the rare gas limit <C 1, we have (n) = 1 + £_/2^ -q = -q + K2rj^, i.e. the increase of the mean 

multiplicity is dominated by the contribution of the second order cumulant K2 — (2roA)~'^. The corresponding 
multiplicity distribution then becomes somewhat wider than the Poissonian one (see Eqs. (|27|), (^)): w(n) — 
const - (l + C^X2)??"/n!. 

In Fig. 2c we demonstrate the approach of the mean multiplicity (ri) to the limiting scaling value = ^/(l — 
though only for ^ very close to the explosion point ^ = 1 > 0.99). Instead, in the region of ^ < 0.9 indicated by 
present experiments, we can see, in agreement with Eq. ( pq ) an approximate ^-scaling of the ratio {n)/rj (Fig. 2b). 

Since, in the realistic event generators, the multiboson effects can be accounted for only in the lowest order cumulant 
approximation |l8| , it is instructive to study the saturation of these effects with the increasing number Ncum of the 
contributing cimiulants. In Fig. ^ we show iVcum-dependence of the ratio (n) /1] of the BE affected mean multiplicity 
to the initial one for different values of the density parameter ^. For example, at ^ = 0.8 this ratio saturates at 
Ncum ^ 10 (^ 40% increase of (n)). At Ncum = 4, representing a practical limit due to the numerical problems 
p8[ , the effect is underestimated by ~ 25% {{n)/ri « 1.3 instead of 1.4). The situation is more optimistic for lower 
densities. Thus at ^ = 0.5 the effect (^ 15% increase of {n)) practically saturates at Ncum — 4. 

Up to now, we have considered the symmetrization effect on the production of only one type of identical bosons. 
For a system of charged and neutral pions, the total symmetrization weight in the model coincides with a product of 
the separate BE weights w„^, a;„_ and LUno provided the pions are emitted in unpolarized and uncorrelated isospin 
states, the effect of their FSI is negligible and there are practically no restrictions due to energy-momentum and 
isospin constraints. The latter assumption may be reasonable at high energies when the subsystem of interfering 
pions represents a small part of the produced multiparticle system. The initial distribution of pion species is then a 
trinomial one, strongly peaked at = rt_ = hq — n/3: w{n^, 7i_, no) = n!/(3"n+!n_!no!). After the symmetrization 
in the large-n condensate limit (a;„. = nj//?"'), it becomes independent of rii and yields a substantial probability 
(2(n — no + l)/[(n + l)(n + 2)]) of any value of no at a fixed total pion multiplicity n. Particularly, the production 
of so called Centauro (anti-Centauro) events containing mainly charged (neutral) pions then becomes possible. The 
probability of the extreme charge configurations can be enhanced even stronger in case of isospin constraints, for 
example, if pions were produced in isosinglet pairs 1 5|,p|,p5[ . The latter mechanism can be of particular importance in the 
near-threshold multipion production,^ due to the limited total isospin and charge. Unfortunately, the n! enhancement 
of the near-threshold condensate production will be more than compensated by the phase space suppression factor of 
(p„/A)3" ~ (n!)-3/2, .^vherep„ = [2m{y/s -Y,^Tn,)/n]^^'^ is the mean pion momentum near threshold. 



B. Single-particle spectra 



The influence of the BE effect on the single-boson spectrum for a given boson multiplicity n, can be seen from Eqs. 
(p5|)-(|30|), (^) and (^9|). At sufficiently large momenta, when the local density (/n)p remains small even at large n, 
this spectrum is dominated by the contribution /3„P(p), /?„ = na;„_i/a;„, of the initial spectrum. In such a rare gas 
limit, ^„,p = ^„exp(-p2/2A2) < 1, we can write (see Eqs. (p5|), (pi) and (p^): 



Nnip) 



n 6-i(p, p) + n(n 

n[l-(n-l)if2]Gi(p,p) 
nP(p) -I- n(n 



1) G2(p,p) 

n(n- l)G2(p,p) 
l)if2[2^/'P(2^/'p)-P(p)]. 



(67) 



A proposal of an experimental study of the near-threshold multipion system at Serpukhov accelerator was recently discussed 
by V.A. Nikitin. Similar idea was also communicated to us by L.L. Nemenov. 
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Otherwise, at large local densities, iV„(p) is determined by the asymptotic large-density spectrum |15|: 

Nn{p) -> n (-^)'^'exp (-^p') = nPc(p), (68) 

associated with the BE condensate and corresponding to the asymptotic (large-n) value ro/(4A) of the slope pa- 
rameters b^. Note that Pc{p) is normalized to unity and that, at /3 = 1, it coincides with the initial distribution 



It is clear from Eqs. (|54|), (p5| ) and (58) that for small momenta, p < A(roA — 1/2)"-'^/^, the condensate regime in 
Eq. ( |68| ) settles on condition ^„ > e. For larger momenta, we must take into account that the condensate contribution 
vanishes much faster than that of BE gas, thus leading to much stronger condition of the condensate dominance: 

^„ > (2roA)~3/2 exp[(2roA - 1)pV2A2]. (69) 

Similar to the fixed multiplicity case, the inclusive single-boson spectrum at small local densities tends to riP{p) = 
Gi(p,p) and, at large ones, it approaches the asymptotic high-density spectrum (see Eq. (|68|)): 

N{p) ^ (n)Pe(p). (70) 

The transfer of the initial spectrum to the high-density one is demonstrated for the inclusive distribution in Fig. 4a 
and, more clearly, for ^ closer to the explosion point f = 1, in Fig. 4b |^ 

Tailing the large- and small-n behavior of the Gn-functions at nt ~ rpA (see Eqs. (|4|) and (p7|)), we can approximate 
the inclusive single-boson spectrum in Eq. ( |6^ ) as 

N{p) = [l + Ip / 2^/2 + . . . + / n,3/2] ^P(p) + C^Pcip) 

^7Vg(p)+7Ve(p), (71) 

where = {27r)^/^riP{p)/ro'^ = fexp(— p^ / 2A'^). Clearly, for large phase space volumes, {tqA)^ ^ 1, the two terms 
in Eq. (|7^) can be interpreted as contributions of the BE gas and BE condensate respectively. Like in the fixed 
multiplicity case, the condensate dominates on condition (|6^) with the substitution n — > (n) (^„ {n)/(3). 

In Fig. 4 we compare the inclusive single-boson spectra with the approximate formula (|7l|). A good agreement is 
obtained despite the calculations were done for not very large phase space volumes. Some underestimation of N{p) 
at intermediate local densities ^p (Fig. 4a) and the corresponding underestimation of (n) at moderate densities ^ 
(Fig. 2b) become weaker for larger systems (larger multiplicities in Fig. 2b) due to increasing number rit ^ tqA of 
the terms in Eqs. ( |66|), (f n\) and thus - decreasing relative contribution of the tailing region. We may conclude that 
the accuracy of Eqs. ( |6q ) and ( |7l| ) is reasonable for the systems produced in heavy ion collisions at SPS and that it 
will be even better for larger systems at RHIC and LHC. Experimentally the effect of BE "condensate" was searched 
for at SPS CERN as a low-pt enhancement, however, with rather uncertain results (see, e.g., p^]). 

It follows from Eq. ( [7l|) that for sufficiently large (nj = tqA ^ 1) and not very dense (^ ^ 1) systems, similarly 
to the ^-scaling of {n)/ri (see Fig. 2b), the ratio N{p) /[i]P{p)] scales with the local density parameter ^p. It appears 
that analogical scaling takes place also at fixed multiplicity n.PI In Fig. 5 we show the ratio of the single-particle 
spectrum at fixed n to the dominant large-p contribution /3„_P(p) of the initial spectrum calculated at p = as a 
function of ^„ for various multiplicities n. An approximate ^„-scaling is seen up to ^„ of the order of unity. At larger 
^„ this ratio approaches the condensate limit (2roA)^/2^„ which no more scales with ^„ (see the corresponding curves 
in Fig. ||b). What scales at large ^„ is not the ratio of the two contributions but the ratio of their integrals, n//3„, 
the limiting value of which is just equal to S,n — n/P since, according to Eq. (p5|), /3„ = nujn-i/^n — > P for n > e/3. 

Note that, in the absence of simple analytical approximations for the fixed-n spectra, the approximate or Cra,p~ 
scaling can be used to overcome technical problems with factorially large numbers at high multiplicities. For example, 
Nn{0) at a given density ^n, < 1 can be obtained by calculating Nn'{0) at a smaller multiplicity n' < n keeping the 
same density = ^„ and then rescahng to A'^^(O) = [Pn/ Pn')Nn'{0), where Pn/Pn' = 1 for n' > ep. 



^These results agree with those obtained in refs. (see also [plJlql and references therein) except for an incorrect conclusion 
that the width of the narrow peak due to the BE "condensate" is of 1/ro. 
"in this case, due to the explicit dependence of the particle spectra on the comphcated BE weights uin, there is no analytical 
approximation similar to Eq. (XTw. 
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C. Correlation functions 



It follows from Eqs. ( p5| ) and (^)-(|59|) that, for a given multiplicity ri, the eorrelation function intercept i?n(p, p) = 
i?„(0) decreases and the correlation function width increases with the increasing n or decreasing momentum p, both 
corresponding to the increasing local density parameter S.n,p- 

In fact, for large local densities (see Eq. (|6^)), the condensate behavior is achieved (see the curves in Fig. 5b) and 
the correlation function tends to unity not only at large but also at small q'^. Indeed, in this limit the normalization 
constant c„ — > 1 (see Eqs. (^8|) and (^5|)) and the nominator and denominator of the correlation function ( p5 [) consist 
of about the same number ^ n? of the condensate terms (ro/TrA)'^ exp [-(4p2 + q2)ro/2A] (see Eqs. (g), ^ and 
(p8|))P] The well known [^-^ lowering and widening of the correlation function with the increasing multiplicity is 
demonstrated, in the considered model, in Fig. 6. 

Note however that in the rare gas limit, ^„^p ^ 1, the change of the form of the correlation function with the 
increasing density is rather weak. Thus, writing in this limit 



A^P(Pi,P2) -n(n-l)^^[Gi(pi,pi)Gi(p2,P2) + Gi(pi,p2)Gi(p2,Pi)] 
+n{n - l)(n - 2)^^-^ [Gi(pi, pi)G2(p2, P2) + Gi(pi, p2)G2(p2, Pi) 

'G2(pi,Pl)Gi(p2,P2) + G2(pi,P2)Gi(p2,Pl)] 



nln 



1)[1 - (2n - 3)X2]P(pi)P(P2) [1 + exp(-fo'q')] 



i{n - l){n - 2)23/2X2{P(pi)P(2i/2p2) 



l + expf-|foV ' 



+P(p2)P(2i/2p,) [1 + exp (-fro^q^ - ^)] }, (72) 

ra"^ = ro'^[l — (2roA)~2], we get for the correlation function intercept an n-independent value close to 2: i?„(0) = 
2(1 — ei). Here we have introduced the density parameter 

= j ■ 2^'^K2 exp(-p2 / 2A2) = 2-^'''£,,^^. (73) 

It should be noted that the correlation function i?„(q) becomes less than unity at intermediate q- values and 
approaches the limiting value of 1 from below. This behavior and also the related suppression of the intercept value is 
caused by the BE correlation effect on the single-particle spectra entering the denominator of the correlation function. 
Sometimes this distortion is corrected for by a special iterative procedure. Its result can be described by a simple 
low-q'^ correction factor p| 

i?r(Pi, P2) = (^n - P„q2)i?„(pi, P2). (74) 

In the low density limit of our model, we have 

A„ = l + £„_!, P„ = e„_i[l-(p/A)2cos2V]/8A2, (75) 

where ip is the angle between the vectors p and q. At small q^ , the corrected correlation function (properly normalized 
to unity at large q^) is then 

i?r(Pi, P2) = 1 + exp(-fo'q') + e„-2 [l + exp (-ffo'q')] (76) 



Of course, the absence of the correlation in the condensate limit at fixed multiplicity has nothing to do with the coherence 
effect which is absent in the considered model. See the footnote after Eq. ( p6| ) and the discussion of the inclusive correlation 
function which appears to be different from 1 at whatever high densities. 

^^The iterative correction procedure is usually used for small-acceptance detectors triggered by the requirement of at least two 
identical pious in the detector. The mixed reference sample then differs from the product of the single-particle spectra, being 
much more influenced by the residual BE correlations. The residual correlations can substantially affect also single-particle 
spectra in the case of a small effective emission volume, e.g., in e'''e~-collisions. There are also other reasons for the low-g'^ 
correction factor, like energy-momentum constraints or presence of dynamical correlations (e^., in jets) which are destroyed 
in the mixed reference sample. For this reason the correction factor similar to that in Eq. ( |74| ) is often introduced as a pure 
phenomenological one with An = 1 and _B„ treated as a free parameter. 
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and can be represented in the usual single-Gaussian form: 



K°''{pi,P2) = 1 + A„ cxp(-r„2q2), 

(77) 

A„ = 1 + 2e„^2, rn'^ = tq^ (l - jen-2) ■ 

We see that with the increasing n the effective interferometry parameters A„ and r„ respectively increase and decrease 
slightly, starting from the zero density values of 1 and tq. 

In the low-density limit, simple Eqs. ([T^), ( [77| ) or directly Eqs. ( |67| ) and ( [72[ ) allow one to determine the radius 
parameter rg by fitting the correlation functions i?^"*" or i?„. At higher densities, however, there is no simple analytical 
expression for the correlation function i?„ and the eventual fit would require the use of rather complicated Eqs. (j5^). 
Another possibility is still a simple single-Gaussian fit at sufficiently small q, giving the effective interferometry 
parameters < 1 and r'^ < rg, both vanishing with the increasing local density. The low-density radius tq and the 
density can then be determined comparing and with the model predictions as functions of tq and A. 

In Fig. 7 we show the intercept as a function of the multiplicity n and the local density parameter ^„^p for 
several values of the mean momentum: p — 0, 0.1, 0.2 and 0.4 GeV/c. As expected, the intercept is practically 
constant at low local densities (small n or high p). As condensate develops, the intercept sharply falls down. The 
sharpness of this drop is however less pronounced at higher momenta even if plotted as a function of the local density 
parameter ^n,p- Clearly, this lack of density scaling is related to a strong decrease of the condensate contribution 
with the increasing momentum. In fact, for the momenta p > A(roA — 1/2)"-'^/^ the low-density parameter ^„_p 
strongly overestimates the local density in the region of the condensate dominance (see Eq. (|69|)). To demonstrate 
the possibility of the observation of the condensate effect for a speckle of a large number of soft pions not following the 
ordinary proportionality rule between the freeze-out phase space volume and pion multiplicity, in Fig. 7 we indicate 
by the arrows the intercept values corresponding to ^„ = 3^ = 1.5 (n « 3(n)). 

In the inclusive case corresponding to the initial Poissonian multiplicity distribution, the correlation function in- 
tercept is equal to 2 for any local densities (see Eq. (^)). At very large local densities the two-boson spectrum 
approaches twice the product of the single-boson ones (see Eqs. (|5^ ) and (|6^)) so that the inclusive correlation 
function tends to the limiting value of 2 even at rather large relative momenta. The corresponding increase of the 
width of the correlation function with the increasing density parameter ^ is demonstrated in Fig. ^. 

Note that at high local densities both the nominator and denominator of the correlation function at small are 
dominated, like in the case of a fixed multiplicity, by the universal condensate terms, their numbers being about 
{n{n — 1)) and (n)^ respectively. The difference between the inclusive and fixed multiplicity correlation functions, 
R ^ 2 and _R„ ^ 1, is due to the fact that at high densities the initially Poissonian multiplicity distribution approaches 
a much wider BE one, for which {n{n — 1)) — 2(n)^. 

Tailing the approximate equations (^) and (^^ for the G„-functions at nt ^ tqA, we can analytically follow the be- 
havior of the function G'(pi, P2) (determining, according to Eq. (|6^), the inclusive two-boson spectrum A^(pi, P2)) sim- 
ilar to Eq. ([n]), modifying it by the substitutions: ^p~^ ^p"^ exp(— q^rQ/2n) and -Pc(p) ^ -Pc(p) exp(— q^ro/4A). 
Particularly, for large systems ((2roA)^ ^ 1), we have at g ^ 0: 

G(pi, P2) = iVg(p) exp(-q2rgV2) + A^e(p) (78) 

where = ro/(2A) and 

nt nt 

r,' = ro' |p"n-^/2 / J] Ip-n'^/^. (79) 

n—l n—1 

Note that at low local densities (^p <C 1) the effective radius rg coincides with "Tq = ro . With the increasing local 
phase space density it slightly decreases. The maximal reduction factor of 1/^/2 is achieved for large (roA ^ 1) and 
dense (^ — > 1, p — > 0) systems. Considering the limit q ^ and rgA 1, we can neglect the g-dependence of the 
condensate term (rc ^ r g) a nd of the product of the single-particle spectra in the denominator of the correlation 
function and, using Eqs. p^ ) and (^, write 



Note that in [g/7| a similar parameterization was used for the uncorrected correlation function. This led to different estimates 
of the interferometry parameters A„ and r„ in the considered model. 
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It follows from Eq. (|8^) that the condensate contribution leads to an additional reduction of the interferometry radius 
squared (defined as a low-q'^ slope of the correlation function) as compared with the case of a pure BE gas. In case of 
a dominant BE condensate the interferometry radius tends to zero whatever large is the geometric size of the system. 

Note however that, due to the non-Gaussian character of the correlation functions at large phase space densities, 
their real width is determined by the large-q^ behavior which reads as: 

i?(pi,P2) = l+( ^XcM^r^W)- (81) 



.^g(p). 

Requiring i?(qcff) = 1 + 1/e, we get for the corresponding effective radius squared: 

reff2 = qeff-2 = re V {1 + 2 ln[A^c(p) / A^g(p)]}. (82) 

Thus, compared with the vanishing of the low-q^ slope of the correlation function (as (1 — ^)), the effective radius 
squared vanishes at large phase space densities much slower (as 1/| ln(l — ^)|) (see Fig. ||). 

It should be noted that Eqs. ( |6^ ) and ( |63| ) assume that the initial Poissonian multiplicity distribution extends to 
any arbitrarily large number of bosons. In reality, however, this number is limited due to the finite available energy 
(for a study of the energy constraint effect on pion multiplicity see second paper in ref . [ p^ ) . It is therefore interesting 
to see how fast the semi-inclusive spectra approach the inclusive limit with the increasing number ri,„ax of the included 
pions. In Fig. 9 we demonstrate the nmax-dependence of the semi- inclusive correlation functions for a fixed value 
of the density parameter ^ = 0.95 and, in Fig. 10 - the rimax-dependence of the correlation function intercepts for 
different ^-values. We can see that the width of the semi-inclusive correlation function increases with the increasing 
f^max, while its intercept decreases at small nmax, reaching a minimum at rimax ~ (n), and then approaches the limiting 
value of 2 roughly as logrimax- The inclusive behavior is practically saturated at a moderate number of the included 
pions riiiiax = k{n), where k increases with the density parameter £ from about 3 at ^ = 0.89 to about 5 at ^ = 0.99.P^ 
Thus the neglect of the energy-momentum constraints in Eqs. ( |62|) and ( |6^ ) can be justified provided (n) ^ ^fsjm, 
e.g., in the usual case of a logarithmic increase of the mean multiplicity with the c.m.s. energy ^/s. The situation 
can change in the case of very large and dense systems dominated by a soft condensate. Then the regime (n) ~ 
can settle, the energy-momentum constraints leading to the reduction of the maximal effective number of produced 
pions to rimax ~ {n) and, as a result, to the suppression of the measured inclusive correlation function. Clearly, such 
an eventual suppression has nothing in common with the coherence effect. 

VI. DISCUSSION AND CONCLUSION 

We have illustrated an approximate scaling of multiboson effects with the density parameters ^,<^ti, • ■ • (see, e.g., 
Figs. I and I). It means that though our numerical results were obtained for typical AGS or SPS multiplicities of 
the order of tens or hundreds of pions, they are approximately valid also for higher multiplicities expected at RHIC 
or LHC energies. 

The value of the density parameter ^ can be estimated with the help of Eq. ( ^ ) relating the phase space density 
in the rare gas limit with the integrated correlation function. Thus using the usual Gaussian parameterization for the 
correlation function in the longitudinally co-moving system (LCMS): 

R{pi,P2) = 1 + Aexp(-r^q2 _ ^2^2 _ rlql), (83) 

where a;, y (y || z x p) and z denote the outward, sideward and longitudinal directions respectively and parameterizing 
the single-particle spectra as 

, . _ dn exp(-(mt-m)/r) 
ay zttJ (1 + ra)mt coshy 



^''The increase of the saturation point nmax/(n) with the density is related with the increasing condensate contribution which, 
for the ideal BE gas, is characterized by very large multiplicity fluctuations. 
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we arrive at the mean pion phase-space density 



/ p —r^N{p) coshy = A— , — r -p, 85) 

V 2 Vi [1 + m)mt ay 



where V = r^ryr^ is the LCMS interference volume.]^ For soft pions {pt ~ and y ~ 0) at SPS energies this quantity 
is typically ~ 0.2. Since this value is sufficiently small, we can compare it with the model phase space density in the 
rare gas limit (/)p=o ~ v/iV^roA)'^ (see Eq. (^7|)) and get ^ « 0.4 — 0.5. For such values of the density parameter our 
calculations point to rather small multiboson effects in the ordinary events. These effects can show up however in the 
events containing sufficiently high density fluctuations. Particularly, the condensate effects could be seen in certain 
high multiplicity events (see, e.g., Fig. [7|) in which the phase space volume or subvolume ~ (roA)'^ does not follow 
the increasing multiplicity (as it presumably does in the ordinary events) and remains sufficiently small to guarantee 
a nonvanishing factor determining the condensate size (see Eqs. (!§) and (|7l|)). 

Since at present energies the LCMS interference volume V seems to scale with dn/dy, the freeze-out of the pions 



occurs on average at approximately constant phase space density (see Eq. (85)). In the rare gas regime, based on 
the density scaling one can then expect about the same relative size of the multiboson effects also at RHIC and LHC 
energies, up to a slight increase in ^ due to the vanishing of the finite-size corrections with the increasing phase space 
volume of the emitting system. At the same time, the growing phase space volume will lead to suppression of the 
average condensate contribution, determined by the factor 

Considering the multiboson effects in the low (BE gas) and the large (BE condensate) density limits, we have 
obtained simple analytical formulae accounting for the finite size of the phase space volume and allowing to follow the 
dependence of the mean multiplicities, single-boson spectra and two-boson correlation functions on the phase space 
density parameters. In principle, these formulae provide a possibility to identify multiboson effects among others. 
Particularly, the width of the low-pt enhancement due to the BE condensation decreases with the size of the system 

— 1/2 

as Tq ' and this narrowing makes the observation of the effect easier. 

The results of the considered simple model should not be taken, however, too literally since: 

a) due to its static character, the model does not explicitly account for the experimental indications on a constant 
freeze-out phase space density and the related expansion of the emission volume. The qualitative application of our 
model to heavy ion collisions is however possible in the limited phase space regions. For example, the pions with 
a rapidity difference grea ter than about unity have to be considered as originating from different static sources; As 



discussed in Section [IB , the residual slow relative motion, if decoupled from other source characteristics, merely 
leads to a wider momentum dependence of the emission function in Eq. (p2[): ^ -I- A§; 

b) due to a mixture of different production processes {e.g., due to contribution of different impact parameters), the 
real multiplicity distribution and particle spectra will be rather weighted sums of those in Eqs. (p9|)-(pl|) calculated 
with different sets of the parameters 77*,rg, A*. As a result, near the condensate limit, the multiplicity distribution 
can be wider than the BE one and the intercepts of the inclusive and fixed-n correlation functions can differ from the 
respective single-process values of 2 and 1; 

c) when estimating the freeze-out phase space density from the experimental data, the multiboson system is con- 
sidered as a homogeneous medium. However, there can be large local density fluctuations - speckles which can give 
rise to noticeable multiboson effects even at a moderate value of the mean phase space density; 

d) on the other hand, the multiboson effects can be somewhat suppressed due to a possible violation of the 
factorization assumption in Eq. ^ or due to the lack of the reflection symmetry of the emission volume. In latter 
case the functions G„(pi, P2) are no more real; 

e) for identical charged pions, the BE effects are also suppressed due to the Coulomb repulsion. Since this repulsion 
is important only in a weakly populated region of very small relative momenta determined by the pair Bohr radius 
a = 387 fm, the suppression of the global BE weights cun is rather small. For example, for uj2 this suppression, being 
about (aroA^)~^, is usually less than one per mill. The Coulomb distortion of the global multiboson effects is therefore 
negligible in the rare gas limit. Nevertheless, since the Coulomb repulsion destroys the formation of the condensates 
made up from positive and negative pions in the disjoint phase space regions, it can lead to noticeable differences 
between charged and neutral pions in dense systems. Particularly, we can expect a decrease of the charge-to-neutral 
multiplicity ratio with the increasing phase space density. 

Because of large numbers of positive and negative pions produced in heavy ion collisions, one could also raise a 
question about importance of the Coulomb screening effects violating standard two-body treatment of the correlations 



better estimate may require the substitution [0] A — > A^^^ in Eq. (|85|). 
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in the low density limit. There are however arguments showing that the screening will be of minor importance even 
at LHC (see also |2^]). Note that in the scenario with a constant phase-space density the corresponding Debye 
radius 

rz, = [47r(p++p_)eVT]-i/2^ (86) 

where — 1/137 and + p_ is the total density of charged pions in the configuration space, will be also constant, 
up to a weak energy dependence due to the temperature T. Assuming that pions with a rapidity difference greater 
than unity come from spatially disjoint regions of phase-space, we can put p5| 

and obtain td k, 15 fm at ~ 0.1 and T k, 200 MeV/c (td ^ '^/T). Thus at LHC energies we can expect the 

characteristic distances between the pion production points comparable or larger than the screening radius leading 
to a suppression of the usual two-particle Coulomb effects. In fact, two charged pions produced at a distance r* > td 
start to feel their Coulomb field only after some time when the density decreases to a value corresponding to Debye 
radius larger than r* . During this time the vector of the relative distance between the pion emission points increases 

by [0 

Substituting r* by r* -I- Ar* in the argument of the Coulomb wave function, we can see however that the suppression 
of the Coulomb effect can be substantial only in the region of large relative momenta k* > (jnTY/'^ where the 
correlations due to QS and FSI are already neghgible (see, e.g., [ |T9| , p8| ). 

Finally, to make easier the understanding of our results in context of other papers and, for the reader's convenience, 
we compare our results (including those in jl^) with the results of refs. Pj6[ p^Jl7| | and the recent papers p9|-p4 36-39|] 
which appeared either after the present work was basically completed or during the process of its evaluation. 

In the pioneering papers [^,^, the analytically solvable model discussed in Section IV was introduced. A simple 
algorithm was given allowing to calculate BE modified multiplicity distribution, single- and two-boson spectra in terms 
of the quantities C^''"** = gn/n, G^''"**(pi, P2) = G„(pi,p2). The details of the calculation of the G„-functions and 
of the two-boson spectra were not given in (|]^. They can be found in [p^|-p^ (see also |2^,^). In refs. (2^ this 
technique was extended to the spectra and correlations of three or more pions. It was shown that, in the inclusive 
case, the old formalism, not accounting for the multiboson effects, can be recovered by a redefinition of the Wigncr-like 
density which then becomes more narrow both in momentum and configuration space. 

Regarding the choice of the Wigner-hke density D{p,x), in the original papers [^,|| and in |l^, it was slightly 
different from that in Eq. (|4^), corresponding to the substitution exp(— q^/2A^) — > exp(— po/T"), where T = A^/m. 
This choice is, in fact, in contradiction with the uncertainty principle, the latter requiring an energy independent 
density in the considered case of a fixed emission time t = Q (see Eq. (po|)). (For the same reason, the freedom 
in the choice of the parameters tq, A is limited by the inequality tqA > 1/2, the equality corresponding to the 
zero distance between the emitter centers.) The incorrect choice of the Wigncr-likc density in [p|,|6|jl^] can be cured 
by the substitution Tq^Tq — (2A)^^. To recover the recurrence relations in Eqs. (|0)jbesides this substitution, 
one has to take into account that the width and the slope parameters used e.g. in |6| are related to ours as: 
i?o = 2rg,Po = o.n = + b~ , gn = 2(6+ — b~) and that a factor of I&js missing in Eq. (66) for gn+i in 



The analytical solution for the G„-functions has been found in ]17| , |30| , |31[ . In |17| the same form of the Wigner-like 
density as in Eq. (^) was used with the parameter s R^ ff = r§ and 2A^. However, since the analytical solution 

was derived based on the recurrence relations of ref. |l6|| , it has also to be cured by the substitution r§ ^ Tq — (2A)~^. 
Particularly, after this substitution, Eq. (200) in the preprint version of for the critical multiplicity rjc (coinciding 
with Eq. (9) in corrected for the misprints ||l^) then reduces to the simple result p^: r/c = (i. In the 
analytical solution given in Eqs. (2.14) corresponds to the emission function in Eq. ([42|) with the parameters 
Tq = i?^/2 + (2po)~^, A^ = Requiring a matrix algebra, this solution is however less transparent compared 



Eqs. (2.14b) and (2.14d) have to be corrected for the misprints by the substitutions; 1/4 (l+c/2)/4 and (1+c) (l+c/2) 
respectively. Note that the other solution given in Eqs. (4.3), (4.4) of ref. presumably contains an error since, corresponding 
to a Gaussian factorizable model, it does not reduce to the general solution in Eqs. ( ^9[ - pl| ) (particularly, it shows no condensate 
behavior at high densities). 
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with that in |T^. In |3^, the analytical solution was obtained using the technique of the density matrix in the 
1-dimensional momentum space. Generalizing this solution to three dimensions (by substituting the normalization 
factors Aq/[1 - (1 - Aq)"] and (27rA^)^^/^ by their cubes) and relating the notation of ref. Hlj] to ours; A^^ = 8&+, 
ip, = 26^, uj^^ = (1 — Ao)" = /o", L{pi,p2) — Gn{pi,P2), one can see that it then coincides with the solution in Eqs. 
&■ 

As for the analytical approximations of the mean multiplicity, single- and two-boson inclusive spectra, similar results 
as in ref. [T5| ] have been obtained in |17| in both the low and high density regimes. The behavior of the single-boson 
spectrum in the low density limit was also studied numerically in - a linear increase of the relative correction with 
the phase space density was found (in agreement with Eq. (71) at small ^p). The condensate behavior at very high 
densities was also obtained in js^. A comment requires the effective radius squared, estimated as ||3l|] r^^j — (2(q^))~^, 
where the averaging over the correlation term is assumed. With the increasing mean multiplicity, r^g decreases from 
Tq to rg/(2roA). Note, however, that at high densities, r^g has little to do with the interferometry radius squared. 
The latter accounts also for the change of the single-boson spectrum and, with the increasing density, decreases from 

- (2A)-2 to 0. 

A special comment requires the normalization of the correlation functions. There are two popular definitions of the 
two-particle correlation function, both representing a ratio of the two-particle spectrum to the product of the single- 
particle ones, but differing in the normalization. In case (I) the spectra are normalized to 1 while in case (II) - to the 
numbers of single- and two-particle counts. This definitions correspond to Eq. (BS) or (Bq) with (I) c„ = n/(n — 1) or 



c = in 



/{n(n — 1)) and (II) c„ = c = 1. In the inclusive case, based on the thermal-type models (see the comment 
after Eq. (|63|)), the second choice was advocated |^,^. Recall that the first choice (I) would be preferable in case 
of a fixed pion multiplicity provided a large phase space volume (see discussion after Eq. (p9|)). Generally, however, 
there is no well-defined normalization and, for a reliable comparison with the experimental correlation functions, a 
free normalization parameter, depending on the production mechanism (not necessarily a single thermal one) and 
experimental conditions, has to be introduced p5[ |. 

The role of the normalization was misunderstood in jl^Jl^, where the decrease of the intercept of the inclusive 
correlation function (I) with the increasing phase space density, obtained in the factorizable Gaussian model, was 
incorrectly interpreted as a coherent laser behavior. Recall that pions in a coherent state would have a narrow 
Poisson multiplicity distribution, while the BE condensate is characterized by very wide multiplicity fluctuations. 
Clearly, the coherent pion production requires a special mechanism (not present in the model), like that leading to 
possible formation of the disoriented chiral condensate - DCC (for a review, see [Q). Possibilities of experimental 
investigations of BE condensate and DCC phenomena have been recentlydiscussed in . A discussion of statistical 
physics aspects of the multiboson effects can be found in [|jl|,|^,||,|3^] . 
In conclusion we summarize the results. 

- The influence of the multiboson effects on boson multiplicities, single-boson spectra and two-boson correlations, 
including an approximate scaling behavior of some of their characteristics with the phase space density {e.g., {n)/ri 
vs ^ or Nn{0) / Pn vs ^„), has been demonstrated using the analytically solvable Gaussian model. 

- The approximate analytical formulae are given allowing to follow the dependence of these quantities on the phase 
space density parameters thus providing a possibility for the identification of the multiboson effects among others. 

- The meaning and the applicability conditions of the model factorization assumption are clarified using the phys- 
ically transparent Kopylov-Podgoretsky ansatz of classical one-particle sources. For heavy ion collisions, the factor- 
izatiom assumption is expected to be valid in case of an impact parameter selection. 

- The lowest order cumulant approximation, suggested for a practical account of multiboson effects in realistic 
transport code simulations ||l8|| , has been shown to be reasonable at moderate densities indicated by the experimental 
data. 

- At high densities, the spectra are mainly determined by the universal condensate term Pc(p) (e-5-, -^^^-'(p) 

(2) 

nPc{p) and Nn (pi, P2) ^ '^('^— l)^c(Pi)^c(P2)) and the initially narrow Poisson multiplicity distribution approaches 
a wide BE one: (n(n — 1)) — + 2(rt)^. As a result, the intercepts of the inclusive and fixed-rt correlation functions 
(properly normalized to 1 at large l^j) approach 2 and 1, respectively and their low-q^ slopes rapidly vanish with 
increasing density; the corresponding increase of the apparent correlation function width is however rather slow - 
logarithmic in the density. 

- It is found that, even near the condensate regime, the inclusive characteristics saturate at rather moderate 
multiplicities of some multiples of (n) thus justifying the neglect of energy-momentum constraints in the considered 
analytical model. The latter are however important near a multipion threshold, particularly making impossible the 
production of a very cold BE condensate. 

- Though spectacular multiboson effects are hardly to be expected in typical events of heavy ion collisions in present 
and perhaps also in future heavy ion experiments, they can clearly show up in certain classes of events containing 
sufficiently high density fluctuations. 
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FIG. 1. The muhiplicity distribution of neutral pions for a) A = 0.25 GeV/c, m = 2.1 fm, = 30 and b) A = 0.25 GeV/c, 
ro = 1.5 fm, ri = 10, where 77 is the mean multiplicity of the initial Poissonian distributions (dotted curves). The dashed line 
corresponds to Eq. (|6l|). 
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FIG. 2. The slope parameter b of the exponential tail c ■ exp( — 6 ■ n) of the multiplicity distribution (a), the ratio of the 
mean multiplicity to the initial Poissonian one (b) and the mean multiplicity (c) as functions of the density parameters ^ = rj/P 
and ~ ?)i ^ ~ 0-25 GeV/c. The curve in (a): 6 = — In^, the curves in (b), (c) are calculated according to the tailing 
approximation in Eq. (^6|), the line in (c): (n) = £,/{! — £,) = u. 
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FIG. 3. The ratio of the mean multiplicity to the initial Poissonian one as a function of the number of the contributing 
cumulants for different values of the density parameter ^. 
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FIG. 4. The inclusive single-particle spectra corresponding to the density parameters a) ^ = 0.89 and b) ^ = 0.95, 0.99 and 
0.998 (the radius ro slightly varies near 2 fm); the corresponding mean multiplicities are a) 64.3 and b) 33.5, 113.7 and 433.8. 
The histograms represent the exact result, the full curves are calculated according to the tailing approximation ^n\i, the dotted 
ones represent the contributions of the two (BE gas and BE condensate) terms in Eq. ( |7l| ) and the dashed curves correspond 
to the rare gas limit r;P(p). 
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FIG. 5. The ratios of the BE affected single-particle spectrum to the dominant large-p contribution PnNi{p) 
Pn = nujn-i/i^n, Calculated at p = as functions of the density parameter ^„ = n//3. The curves represent the large-^„ 
limit: (2roA)3/2{„ = [2{n/^„y/^ - if^^n- 
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FIG. 6. The two-pion correlation functions for the multipUcities increasing from n = 2 to 362 with a step of 60 (the 
corresponding density parameter ^„ ranging from 0.04 to 7.2 with a step of 1.2). The higher is the multiplicity the lower is the 
intercept of the correlation function and the larger is its width. 
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FIG. 7. The intercept of the two-pion correlation functions as a function of the multiplicity n and the density parameter 
$n,p for several values p — 0, 0.1, 0.2 and 0.4 GeV/c of the mean momentum of the two pions. The arrows on the interpolating 
curves indicate the intercept values corresponding to = 3^ = 1.5 (n « 3(n)). 
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FIG. 8. The inclusive two-pion correlation functions demonstrating the increase of the correlation width with the increasing 
density parameter ^. The different ^-values are achieved by slight variations of the radius ro around 2 fm. 
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FIG. 9. The semi-inclusive correlation functions including the pion multiplicities from to rimax for different values of rimax. 
The dotted curve is the inclusive (rimax oo) correlation function. The conditions are the same as in Fig. ^ for the density 
parameter ^ — 0.95. 
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FIG. 10. The intercepts of the semi-inclusive correlation function including the pion multiplicities from to rimax as functions 
of Umax for different values of the density parameter ^ — 0.89, 0.95, 0.99 and 0.998; the arrows indicate the corresponding mean 
multiplicities (n) — 22.3, 33.5, 113.7 and 433.8. The conditions are the same as in Fig. ^. 
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